Spinor Dynamics of Quantum Accelerator Modes near Higher Order Resonances. 
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Quantum Accelerator Modes were discovered in experiments with Kicked Cold Atoms in the 
presence of gravity. They were shown to be tightly related to resonances of the Quantum Kicked 
Rotor. In this paper a spinor formalism is developed for the analysis of Modes associated with 
resonances of arbitrary order q > 1. Decoupling of spin variables from orbital ones is achieved by 
means of an ansatz of the Born-Oppenheimer type, that generates q independent band dynamics. 
Each of these is described, in classical terms, by a map, and the stable periodic orbits of this map give 
rise to quantum accelerator modes, which are potentially observable in experiments. The arithmetic 
organization of such periodic orbits is briefly discussed. 

PACS numbers: 05.45.Mt, 03.75.-b, 42.50.Vk 



INTRODUCTION. 



A kicked system is a Hamiltonian system that is periodically driven by pulses of infinitesimal duration. More 
than thirty years after the invention of the paradigmatic model of such systems, namely the Kicked Rotor (KR) 
' Kicked Quantum Dynamics is still the focus of active research for a twofold reason. On the one hand, it has 

given birth to an ever increasing list of variants of the basic original prototypes, which have provided formally 
simple models for the investigation of quantum-classical correspondence and of some general properties of quantum 
transport. These include dynamical localization [H, anomalous diffusion [1, [3, decay from stable phase-s pac e 

islands @, H, @] , electronic conduction in mesoscopic devices [13, EH EH , nondispersive wave packet dynamics [l3l | , 
Q" 1 effects of dissipation on quantum dynamics [13], and lately directed transport [HI, Htj Jl7l - On the other hand, 
renewed interest on the physical side has been stimulated by experimental realizations [l8t [T9L I20I [2~fl ] , which are 
now possible, under excellent control conditions, thanks to the science and technology of cold and ultra-cold atoms. 
. Unexpected advances of the theory have been prompted by such experiments. For instance, the so-called Quantum 
' Accelerator Modes (QAM) were discovered in experiments with cold atoms in periodically pulsed optical lattices 
IO i [22I I23I l24l . I25I ]. Their underlying theoretical model is a variant of the Kicked Rotor model, the difference being 
that, in between kicks, atoms are subject to gravity. When the kicking period is close to a half-integer multiple 
of the Talbot time 26], which is a natural time scale for the system, a fraction of atoms steadily accelerates away 
from the bulk of the atomic cloud, at a rate and in a direction which dep end on various parameter values. Though 
QAMs are a somewhat particular phenomenon, their theory [13, [H, S H3] is a vast repertory of classic items of 
classical and quantum mechanics. QAMs are rooted in subtle aspects of the Bloch theory, and have a relation 
to the Wannier-Stark resonances of Solid-State physics [l|. They are a purely quantal effect, and yet they are 
explained in terms of trajectories of certain classical dynamical systems, by means of a "pseudo-quasi classical" 
approximation, where the role of the Planck constant is played by a parameter e, which measures the detuning of the 
kicking period from a half-integer multiple of the Talbot time. This theory hinges on existence of a "pseudoclassical 
■ limit" for e — > 0. That means, for kicking periods close to half-integer multiples of the Talbot time, the quan- 
tum dynamics may formally be obtained from quantization of a classical dynamical system, using e as the Planck 
constant. This system is totally unrelated from the classical system that is obtained in the proper classical limit h — > 0. 



Experimental and theoretical investigation on QAMs are currently focused on novel research lines: the observation 
of QAMs in a Bose-Einstein condensate [3ll . [32l . [HI, which allows a precise control on the initial momentum 
distribution, the analysis of QAMs for special values of the physical parameters [13, HH and, in particular, when 
the kicking period is close to a rational multiple of Talbot time [3H [33] • Theoretical aspects concerning the latter 
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FIG. 1: Momentum distributions, in the time dependent gauge, after t — 100 kicks, for different values of the kicking period 
near the resonance r res = 3n. Red color corresponds to highest probability. The vertical dashed line corresponds to the resonant 
value. White full lines show the theoretical curves (|30[) . with: (left) q — 2, (r,s) — (1, 1) and (right) q = 7, (r, s) = (4, 1) close 
to the higher order resonance r rcs /2ir — p/q — 11/7. The initial quantum distribution is a gaussian wave packet, reproducing 
the experimental conditions. The other parameters are: k = 1 and g = 0.0386. All our numerical simulation refer to the choice 
V{6) = fccos(9. 

problem are considered in the present paper. 

QAMS are connected with an important feature of the KR model, namely, the KR resonances [38j], which occur 
whenever the kicking period is rationally related to the internal frequencies of the free rotor. The dynamics of the 
rotor at a quantum resonance is invariant under momentum translations by multiples of an integer number. The 
least positive integer q, such that translation invariance in momentum space holds, is the "order" of the resonance 
(sect. [TTJ) . The half- Talbot time in atom-optics experiments is the period of the KR resonances of order q = 1 (i.e. 
"principal" resonances), so the originally observed QAMs are related to KR resonances of order 1. 

In this paper we consider quantum motion in the vicinity of a higher order KR resonance (q > 1) in the presence 
of gravity. Numerical (see figll]) and heuristic indications [37f suggest that higher-order KR resonances, too, may 
give rise to QAMs. This has been substantiated by a theory [36[ based on a nontrivial reformulation of the original 
pseudo-classical approximation. It has been remarked that, in the case of higher-order resonances, no pseudoclassical 
limit exists and similarity to the case of quasi-classical analysis for particles with spin was noted, but not explored. 
About the latter general problem [H, H^, it is known that, although no single well-defined classical limit exists, and 
so no global quasi-classical phase-space approximation in terms of a unique classical Hamiltonian flow is possible, local 
quasi-classical approximations are nevertheless still possible, as provided by bundles of trajectories which belong to a 
number of different Hamiltonian systems. 

In this paper we develop a formulation of the problem of QAMs near higher-order resonances in spinor terms. The 
quantum evolution at exact resonance is described by a multi-component wave function, that is by a spinor of rank q, 
given by the order of the resonance [H, Hl| , and is generated by a time-independent spinor Hamiltonian [H, EH • We 
show that the small-e analysis of quantum dynamics is formally equivalent to semiclassical approximation for a particle 
with spin-orbit coupling. Thus QAMs near higher order resonances constitute a particular, though experimentally 
relevant, model system, in which this crucial theoretical issue can be explored. The semiclassical theory in (40j is not 
directly applicable here, because the dynamics is not specified by a self adjoint spinor Hamiltonian, but by a spinor 
unitary propagator instead. We therefore resort to an "adiabatic" ansatz, which allows decoupling spin dynamics 
from orbital motion. In this way we obtain q distinct and independent orbital one-period propagators. Each of them 
may be viewed as the quantization of a formally classical dynamical system, given by a map; however, the "pseudo- 
Planck constant" e explicitly appears in such maps, in a form that precludes existence of a e — * limit for the maps 
themselves, except for the q = 1 case, in which the pseudo-classical theory of ref. [27|, l28| is recovered. 
QAMs, detected by numerical simulations of the exact quantum dynamics near higher order resonances, tightly 
correspond to stable periodic orbits of the maps. The acceleration of the modes is expressed in terms of the winding 
numbers of the corresponding orbits and of the order of the resonance. Moreover, we derive some theoretical results, 
which generalize those obtained in (28l . |29| for the principal resonances: a formula for the special values of quasi- 
momenta, which dominate the mode, and a classification of detectable modes by a Farey tree contruction [3|4^|, as 



a function of the gravity acceleration. 

The paper is organized as follows. In sect|TT]the Floquet operator, describing one-step evolution of a kicked atom 
in a free-falling frame, is recalled and the resonant spinor dynamics in Kicked Particle (KP) model is briefly reviewed; 
in sect IIII1 the quantum motion in the vicinity of a resonance of arbitrary order is related to the problem of a particle 
with spin-orbit coupling. In sect II VI a "formally" classical description of the orbital dynamics, associated to the 
QAMs, is achieved. Finally, in sect fVl connections between the theoretical results and possible experimental findings 
are discussed. 



II. BACKGROUND. 

A. Floquet operator in the "temporal gauge". 

In the laboratory frame, the quantum dynamics of the atoms moving under the joint action of gravity and of the 
kicking potential, is ruled by the time-dependent Hamiltonian (expressed in dimensionless units): 

H L (t) = ^--x + kV(x) V 8(f-nr), (1) 

Z T z — ' 

n— — oo 

where p and x are the momentum and position operators. The potential V(x) is a smooth periodic function of 
spatial period 2ir. Denoting M, T, K,g and 2-k/G the atomic mass, the temporal period of the kicking, the kicking 
strength, the gravity acceleration and the spatial period of the kicks, respectively, the momentum, position and mass 
of the atom in (flj are rescaled in units of HG, G _1 and M. The three dimensionless parameters k,r and r\ in |(TJ|, 
which fully characterize the dynamics, are expressed in terms of phy sical quantities by k — K/7j, rj = MgT/(hG) and 
r = HTG 2 /M = AttT/Tb- T b = 4ttM/(HG 2 ) is the Talbot time ^26] and g = rj/r is the rescaled gravity acceleration. 
Throughout the following h — 1 is understood. 

For j] = 0, the Hamiltonian (fT]) reduces to that of the Kicked Particle model (KP), which is a well-known variant 
of the Kicked Rotor model (KR), corresponding to the particular choice: V(x) = kcos(x). The KP differs from the 
KR because the eigenvalues of particle momentum are continuos while those of the angular momentum of the rotor 
are discrete. Due to Bloch theorem, the invariance of the KP Hamiltonian, under space translations by 2-7T, implies 
conservation of the quasi-momcntum (3, which, in the chosen units, is the fractional part of the momentum. The 
particle momentum is decomposed as p = N + (3 with JV £ Z and < /3 < 1. Conservation of quasi- momentum 
enables a Bloch- Wannier fibration of the particle dynamics: the particle wave function is obtained by a superposition of 
Bloch waves, describing the states of independently evolving kicked rotors with different values of the quasi-momentum 
(called /3-rotors). 

A remarkable feature of Hamiltonian {JJ is that, unless rescaled gravity g = rj/r assumes exceptional commensurate 
values, the linear potential term breaks invariance under 2ir space translations. Such an invariance may be recovered 
by going to a temporal gauge, where momentum is measured w.r.t. free fall. This transformation gets rid of the linear 
term and the new Hamiltonian reads [28| : 

-, +oo 

H a {t) = -{N + p+ r kf + kV{Q) S(t-nr). (2) 

n— — oo 

where 9 — x mod(27r), N = —id/d0 with periodic boundary conditions. 

The quantum motion of a /3-rotor in the "temporal gauge" (that is, "in the falling frame") is described by the 
following Floquet operator on L 2 (T) ( T denotes the 1-torus, parametrized by 9 E [— tt, 7r[): 

Upln) = e -^v(o) e - ti (N+p + nn +i f (3) 
where n € Z denotes the number of kicks. The operator ([3]) describes evolution from time t = tit to time t — (n + 1)t. 

B. Quantum Resonances. 



We consider the problem of Quantum Accelerator Modes in the vicinity of a generic resonance of the /3-rotor. The 
concept of quantum resonance (QR) is reviewed in this subsection. 



A QR occurs whenever quantum evolution commutes with a nontrivial group of momentum translations. A mo- 
mentum translation N — ► N + i (recall H = 1) with I £ Z is described by the operator T e — e l£e . In the following we 
assume 77 = and then the operator §3§ is time-independent. It commutes with T l if and only if [461 ]: i) t/2tt = p/q 
with p, q coprime integers; ii) £ — rq with r£l; Hi) (3 = v jrp + rq/2(mod 1), with v £ Z. In this paper we restrict 
to "primary" resonances, i.e. to resonances with r = 1 and I = q; in this case, g defines the order of the resonance. 
QR of order 1 are called "principal resonances" . 

A theory for QAMs in the vicinity of principal resonances was proposed in [2?l |2S| . In this paper we consider 
quantum resonances of arbitrary order q > 1. The resonant values of the kicking period (i), expressed in physical 
units, coincide with rational multiples of half of the Talbot time. 

We generically denote U roB the operator at resonance, and (3q the resonant values of quasi-momentum, given by 
the above condition (Hi). 

C. Bloch theory and spinors. 

Translation invariance under T q enforces conservation of the Bloch phase £ = 9 mod 2ir/q, taking values in the 
Brillouin zone B = [—n/q, n/q[. Loosely speaking, this means that 8 only changes by multiples of 2ir/q, so £ has the 
meaning of "quasi-position" . As we show below, a Bloch- Wannier fibration of the rotor dynamics holds with respect 
to the quasi-position £, at a QR. 

We use a rescaled quasi-position •& = g£, and accordingly resize the Brillouin zone to [— 7T, tt[. In all representations 
where quasi-position is diagonal, the state \ip) of the rotor is described by a q-spinor 0, specified by q complex functions 
01 W = (^)'IV')) = 1 , ■ ■ ■ , q) ■ We shall use a representation where the spinor which corresponds to a given 

rotor wavefunction ^(0) = (^IV')i is defined by: 

M#) = 4= V^ + m^e™ 19 Z = 0,...,g-1. (4) 

where V'(n), n £ Z, are the Fourier coefficients of ^(0). Equation ([?]) defines a unitary map a of L 2 (T) onto L 2 (T) (giC 9 . 
Under this map, the (angular) momentum operator N is tranformed to: 

N = -id e -> a(7V)a" 1 = -iqdo ® i + J <g> S, (5) 

where 2 and I are the identity operators in L 2 (T) and in C q respectively, and S is the spin operator in C q : 

9-1 

S = X)WI- (6) 

1=0 

where \l), I = 1, ... ,q — 1 is the canonical basis in C 9 . Thus, in spinor representation, the momentum operator is 
the sum of the orbital operator —iqd$ I and the spin operator I ® S. In this picture, the rotor is characterized by 
"orbital" observables ($, —id$) and by the spin observable. 
Bold symbols denote vectors in C 9 and matrices. 

D. Resonant spin dynamics. 

At resonance, quasi-position is conserved under the discrete-time evolution defined by ([3|), so, whenever it has a 
definite value no "orbital" motion occurs, and spin alone changes in time. Therefore the evolution is described by 
a unitary q x q matrix A(z9) such that, as ip(0) evolves into U res tp(0), the corresponding spinor <^(i9) evolves into the 
spinor A(i9)^>(i9). The explicit form of the spin propagator A(#) is easily computed by using ^ under resonance 
conditions. With the specific choice V(9) — cos(0), one finds (details can be found in appendix [A]) : 



A(0) = e^e"* 6 , (7) 
G == G. M = ^(S + poif, (8) 

V(tf) - \ |E C10<i + 1| + |i + l>(i|) + |0)(g - l|e w + |g - l>(0| e ^| . (9) 
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FIG. 2: Eigenvalues of the resonant Hamiltonian H rcs (i?) for different values of the kicking constant k — 1 (red), 3 (green) and 
5 (blue) and (a) q = 2,p — 3, (b) q = 7,p — 11. 

E. Bands. 

The "resonant Hamiltonian" H rcs ($) is a hermitean matrix of rank q such that: 

A(tf) = e - iftres ^. (10) 

It is uniquely defined, under the condition that its eigenvalues (i.e., the eigenphases of A(i9)) lie in [0, 2ir[. Explicit 
calculation of eigenvalues and eigenvectors of A($), hence of the resonant Hamiltonian, is trivial for q — 1, and is 
easily performed for q = 2 in terms of Pauli matrices (43j (such a case is reviewed in appendix [Bj. However, for q > 2 
analytical calculation is prohibitive. 

Eigenphases of A(-#) are smooth periodic function of quasi-position As varies in [— ir, tt[, they sweep bands in 
the quasi-energy spectrum of the resonant evolution described by U rBS [47| . They also depend on the kicking strength 
k and will be denoted by uoi — cu^tf, k) in the following (I = 0, q — 1). In the case q = 1, wo($, k) = fccos(i9). For 
q > 1 the eigenvalues are nontrivial functions of the kick strength k. For fixed q > 2 bandwidths tend to increase with 
fc, eventually giving rise to complex patterns of avoided crossings. Examples of i?- and fc-dependence of eigenphases 
are shown in figJ2]for q = 2 (a) and q ~ 7 (b). For q > 2 the bandwidths depend also on I. 

In the resonant representation (i.e. in the representation in which the resonant propagators (|10p are diagonal), the 
spinor components ((!]) evolve independently. 

III. NEAR-RESONANT DYNAMICS AND SPIN-ORBITAL DECOUPLING. 

We are interested in quantum motion, described by Q, in the vicinity of a QR, namely when the kicking period is 
t = 2np/ q + e, where the detuning e of the period from the resonant period r rcs = 2i:p/q is assumed to be small. The 
one-step evolution operator ([3]) may be factorized as: 

U {n) = U ICS ■ U nr {n) (11) 
[> nr (n) = e-*[^ 2 +^], (12) 

where D n = t(/3 + r\n + r\j '2) — 2-Kpf3o/q. 

A. "Adiabatic" decoupling of spinor and orbital motions. 

Translation invariance (in momentum) is now broken by U nr , so quasi-position is not conserved any more. The 
evolution of a spinor 4> £ L 2 (T) <g) C q is ruled by the time-dependent Schrodinger equation: 

r) +°° 
ie-cj> = H(t)4> , H(t) = eH ICS - £ S(t - n) + H (t), (13) 

n— — og 

H (t) = ^e 2 (-iqdi>®i + i®S) 2 + eD [t] (-iqd#®i + i®S). (14) 



where [t] denotes the integer part oft and (_ff rcs 0) = H res (i?)^(i?). Note that Ho(t) is constant in between kicks. 

Both sides of the Schrodinger equation have been multiplied by e, to make it apparent that the detuning e plays 
the role of an effective Planck constant in what concerns the motion between the <5-kicks. 

The spinor components are mixed by Q14p during the evolution; we use an ansatz of Born-Oppcnhcimcr type in 
order to decouple orbital (slow) motion from spin (fast) motion. 

The detuning e controls as well the separation between the different time scales of the system. At exact resonance 
(i.e. e = 0) the decoupling is exact, because motion is restricted to the eigenspaces of the resonant propagator (fTU)) . 
These subspaces are defined by the spectral decomposition of the resonant Hamiltonian, which we write in the form: 

9-1 

H res W = £ Wi (0,fc)Pi(0) > Pj(0) = \<PM){<PM\, (15) 
j=o 

where <pj ($) is the normalized eigenvector of H rcs (i?) which corresponds to the eigenvalue u>j (•$, k) . For each value of 
the operators Pj($) in (fT5|) are projectors in C 9 . We denote by Pj the projectors in the full Hilbert space L 2 (T) (g>C 9 , 
which act on spinors according to (Pj <£)(•&) — Pj The subspaces Hj whereupon the Pj project are the "band 

subspaces" and are not invariant for the full Hamiltonian (|13p . By using the ansatz that "band subspaces" be almost 
invariant for small e, we next decouple the (assumedly "fast") spin variables from the orbital ("slow") ones. We 
assume that the decoupled evolution inside the band subspaces provides a good description of the exact evolution 
when e is small, because the leading error terms are linear in e. 

Our approximation consists in replacing the exact dynamics, ruled by the Hamiltonian in (|13p . ()14ll by an "adiabatic" 
evolution, generated by the Hamiltonian: 

q— 1 +oo q—l 

H—(t) = Pj H(t) Pj = etf ros • L P i Ho® Pj • (16) 

j — n— — oo j— 

In the case of time-independent Hamiltonians, such projection on "band subspaces" , aimed at separating fast and 
slow time scales, is basically a Born-Oppenheimer approximation [48| . In the case of kicked dynamics this projection 
should be performed on the "effective" , time-independent Hamiltonian iJ off , which generates over a unit time the 
same evolution as does the kicked Hamiltonian. The effective Hamiltonian is not known in closed form, although it 
can be expressed by a sum of infinite terms, ordered in powers of e [43l [49j. Our ansatz is somehow related to a 
rough approximation H ott ~ eZf lcs + Hq- We assume this is valid in some restricted parameter regimes (see further 
comments in Sect. II VI) . 

A spinor m Hj has the form V(#)y>j(i?) with ip G L 2 (T) and may thus be described by a scalar wavefunction ip('d) 
(the amplitude of the spinor on the j-th resonant eigenstate). Evolution inside the band subspace Hj is ruled by the 
"band Hamiltonian" Hj(t) = Pj H(t) Pj and direct calculation by using (113|) shows that band Hamiltonians have the 
following form: 

+oo 

Hj(t) = eu; 3 (d,k) *(*-*') + H^(t), (17) 

t' — — OO 

where: 

H$\t) - -\e 2 q 2 dl - (e 2 q 2 {<Pj\<Pj) + i<?qSj + ieqD [t] ) d„ + \<? (s? - q 2 (<Pj\((>j) - i2qS{) + 

+eD [t] (Sj-iq{ipj\<Pj)), (18) 
where dots denote derivatives with respect to i?, and 

SM = (^WISI^-W) , Sj (tf) = ( Vj (tf)|S|^ W), <W = (v»iWI§ a |^-W) • (19) 

B. Band Hamiltonians. 



We now note that the problem can be formulated as the evolution of a particle in a fictitious magnetic field, which 
takes into account the average effects of spin degree of freedom on the orbital motion. 



We derive a simpler form for the band Hamiltonians feas.(|24p). By the introduction of magnetic vector and scalar 
potentials, the operator (|18|) may be written in the form: 

Hj(t) = \<i 2 q 2 {-id* - A0)f + eqD [t] (-id* - + \^B0) . (20) 

The "geometric" vector potential Aj{d) and the scalar potential Bj(i!)) are determined by the structure of the resonant 
eigenvectors via the following relations: 

A j (#)=i{< Pj w\<p j m-~ s M ( 21 ) 

Bjtf) = S$'(0) + 2q Q^(tf) - q 2 A 2 (0) + ?{4>i{#)\W))- ( 22 ) 

Reality of such potentials follows from IT9|) and from the fact that (<pj('&)\(pj('&)) is purely imaginary thanks to nor- 
malization. The vector potential is gauge-dependent; eigenvectors Pjifi) are determined up to arbitrary ^-dependent 
phase factors and so operator (j2"0")l may be further simplified by a gauge transformation, (pj($) — > (pj{d)e lXj ^\ Under 
such a transformation, Aj{d) changes to Aj{$) = Aj("&) — Aj(i?), and Bj{"d) does not change. The transformation 
may be chosen so that 

Aji'ff) = const. = — 7j i9 — <; = aj , (23) 

where: 

This immediately follows from ([21) and from the requirement, that eigenvectors be single-valued. Note that ^Jj, q is 
the geometric (Berry's) phase [13, [H|. We thus assume A,- = ay: this choice corresponds to the Coulomb gauge. 
In conclusion, in the j-th band subspace, the band dynamics is described by the following Schrodinger equation: 

oo 

Hj(t) = eujj^k) Kt-n) + -e 2 q 2 {-id i) -a, J ) 2 + -e 2 B0)+eqD [t] {-id § ~ a] ) (24) 

71— — OO 

The multicomponent Schrodinger equation (|13p . for the g-spinor wave function (/>(■$, i), is then reduced to g scalar 
Schrodinger equations (j2"4"]l , each of which determines the independent evolution of a rotor wave function ?/>($,<:). 

IV. PSEUDO-CLASSICAL DESCRIPTION OF ORBITAL MOTION. 

We now derive a description of the dynamics of the orbital observables —ids), restricted inside each of the band 
subspaces Hj , by "formally" classical equations of motion. 

We introduce a "pseudo-classical" momentum operator /, defined as follows: 

I = -iedn , (25) 

which differs from the orbital momentum because of the replacement of the Planck constant (= 1) by e. If the same 
role is granted to e in eqs. ([2"4")) . then, in classical terms, the effective band dynamics in the j-th band subspace looks 
like a rotor dynamics, with angle coordinate $ and conjugate momentum /, ruled by the kicked Hamiltonian: 

t— — oo 

Fj{0,I,t) = \q 2 I 2 +D m I-eq 2 a j I+\e 2 B0). (26) 

Terms, independent of / and have been neglected. This Hamiltonian describes a classical kicked dynamics. By 
dropping terms beyond first order in e, the map from immediately after the n-th kick to immediately after the (n+l)-th 
kick is: 

#n+i = $n + q 2 I n + 2-Kflqn + g , mod (2n), 

I n+ i = J„ - eu)j(§ n+1 ,k) , (27) 



'i'.'I- p.H : ':a : 'M'. ■■■] : c : r- -fe^K dU— f=N'i i 1 I — i. ti 
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FIG. 3: Phase portraits of maps (J35J) with k = 1 and g = 0.0386 and (a) t/2tt = 1.445 (e = k = -0.2827), r?? = 2ttS1 = 
3.2261, (r,s) = (1, 1), of map fl33) with j = 1; (b) t/2tt = 1.5025 (e = k = 0.0157), tt? = 3.4401, (r,s) = (23,21), of map 
(Jill with j = 0; (c) p/g = ll/7,r/27r = 1.5375 (e = k = 0.2132), rr? = 3.6023, (r,s) = (4,1), of map with j = 3; (d) 
p/g = 22/15, t/2tt = 1.485 (e = k = 0.1152), rrj = 3.3605, r/s = (8, 1), of map ([28j with j = 15. The periodic orbit has period 
r = 132 and it is associated to 132 stability islands; r and s are not coprime (r — 1056, s = 132). In the inset a magnification 
of one of the small islands of the chain is shown. 

where — rjr / '(27r) and g = q(—eqaj + 7rf2 + r/3 — 2Ttpf3o/q). 

The meaning of the pseudo-classical map (|27|) as a description of the nearly resonant quantum dynamics will be 
discussed in section ITV CI 

A. Pseudo-classical maps and Quantum Accelerator Modes. 

We now describe how quantum accelerator modes appear in the present framework. 

The explicit dependence on time of map (|27p is removed by changing the momentum variable to J n = q 2 I n + 
2TtQqn + g. In the variables (J, $) the map is 27r-periodic in J and so it may be written as a map on the 2-torus: 

i? n+ i = i9 n + J„ mod (2tt), 

J«+r = Jn - £9 2 Wj(l?, fc) + 27r^q mod (2tt). (28) 

In case q = 1, this map reduces to one which was introduced in (2?j], in order to explain the QAMs that had been 
experimentally observed near principal resonances. For q > 1, it has q different versions, labeled by I = 1, . . . , q. Like 
in the case 5=1, the stable periodic orbits of each of these versions are expected to give rise to QAMs. Indeed, each 
stable periodic orbit of map (|28|) corresponds to a stable accelerating orbit of map ([27]) , because the difference between 
momentum /„ and momentum J n linearly increases with time. More precisely, let (i?o, Jo) be initial conditions for a 
periodic orbit of period s and winding number r/s. The increment of J after time ns (measured in the number of 
kicks) is 27rrn; therefore, the increment of the original momentum variable is: 

2tt ( r \ 

Isn - Io = aisn , a/ = — A , (29) 

q \qs J 

with Jo = (Jo — Q)/q 2 - This formula (PTj| yields the acceleration of a stable orbit of the pseudoclassical dynamics 
(|27|) . and it is precisely such orbits that may give rise to QAMs in the vicinity of resonances of arbitrary order. As 
a matter of fact, numerical simulations reveal QAMs near higher order resonances, in correspondence with periodic 
orbits of maps ((25)) . In sect. [V] we explain how the analysis of the stable periodic orbits of maps l|28p may help to 
resolve the complex pattern of QAMs presented in figH] 



Thanks to © and (I2T)]) , the physical momentum N is related to / by N — ql/e+j; therefore, the physical acceleration 
is given by: 

-t(s-°)- (30) 

Although the analytical derivation of the maps is based on the resonant Hamiltonian, which is known in closed 
form only for q = 1, 2, the practical use of (|28p only requires the resonant eigenvalues, which can be easily computed 
by a numerical diagonalization of a q x q matrix. 

In fig [3] examples of phase space of maps (|28[) are shown for q = 2 (a,b), q — 7 (c) and q = 15 (d), in parameter 
regimes in which QAMs are present. The plotted periodic orbits correspond to some of the modes shown in figlU For 
instance, in figE] (a) the stability island of a fixed point of one of the maps |35|) is plotted for t/2tt — 1.455; this fixed 
point corresponds to the huge mode on the left side of figHJ A distribution of phase-space points, which initially fall 
inside the stability island, describes an ensemble of atoms generating the QAM. 

Classical structures, like stability islands, may affect the quantum system only if their size is comparable with the 
effective Planck constant e. The map (among the q set of eq. (|28p). that crucially contributes in determining the 
observed QAMs, is generally the one with the widest bandwidth. 



B. Special values of quasi- momenta. 

A QAM arises when the initial wave packet is centered in momentum around Nq, related to Iq by 

N Q = q— + j = —(J Q + 27m) + gay - - (nfl + t/3 - 2t™/3 /<7) +3 (31) 
e qe e 

with n G Z. As in the case for main resonances, we expect that the modes will be especially pronounced when 
quasimomentum is fine tuned: in view of pip such optimal values of [3 are determined by the condition: 

e , Jo + 27TTO 77 , „ , „ 

0u = --(No - J ~ qoij + A,) + ~ + A) mod (1), (32) 

r qr Z 

with /3o = ^ + § and v = 0, 1, ..,p — 1. A wave packet initially localized in No + @ u will be mostly captured inside a 
QAM; indeed in this case, the overlap between the stabilit y is land and the initial wave packet is maximal. Formula 
(f3"2"| is a generalization of the result derived for q = 1 in [27| and experimentally verified in [3l|; it reduces to the 
expression in [27[ for o-o = (see appendix IB ip . 

This picture is confirmed by figdl in which the quantum phase-space evolution of a /3-rotor, with a quasi-momentum 
given by (I32p . and the pseudoclassical motion are compared. The initial state of the rotor is a coherent wave packet 
centered in the (r, s) = (1, 1) fixed point, plotted in fig[3Ja), corresponding to the e-classical accelerator mode on the 
left part of figlU in the vicinity of q = 2 resonance. The mode moves with an acceleration equal to 0.2988, according 



C. Validity of the Pseudo-Classical Description. 

We now come back to the meaning of the pseudo classical description, as "pseudoclassical" dynamics (|27l) still 
explicitly retains the "Planck constant" e. In the case when q = 1, there is a single resonant eigenvalue, given 
by Wo(i?, k) = fccos($), so the pseudoclassical dynamics (|2T1) has a well defined limit for e — > 0, k — > oo, fee — > k 
with |fe| < oo . This limit dynamics was discovered and analyzed in [2?], [28j]. This is no longer true when q > 1 
and then the relation between the band dynamics and the "pseudoclassical" dynamics ([2"T]) is less transparent. The 
quantum band-dynamics is still, formally, the quantization of the classical kicked dynamics (|27p using e as the 
Planck constant. Nevertheless, the latter dynamics contains the "Planck constant" e in crucial ways, which preclude 
existence of a limit for e — > 0. To see this, note that k) depends on its arguments only through the real variables 
u = ksm(-d/q),v — kcos($/q) (cf. the form of the resonant evolution (IA3|) in appendix fA]) that is, u>j(i},k) = G(u,v), 
where G is a smooth oscillatory function, independent of k. Hence, 



ewj(#,k) = — {cos(&/ q)d u G - sin(tf/ q)d v G} 



(33) 




FIG. 4: Contour plots at times t = 100 (a) and t — 200 (b) of the Husimi distribution of the wave packet of a /3-rotor 
with ft — 0.1672, given by (137[l with j = 1 and v = 0. The rotor is initially prepared in a coherent state centered in the 
(r,s) — (1, 1) fixed point of figEIa). The black spots in the centers of the contours are an ensemble of classical phase points, 
initially distributed in a circle of area ~ e centered at the mode. They evolve according to the e-classical dynamics (|36[) with 
q = -0.5709. The other parameter values axe k = l,r = 1.455 * 27r(e = -0.2827) and g = 0.0386. 



Existence of a limit demands ek — > k; but then, except in the trivial case k = 0, the arguments of the G functions in 
(|33f diverge and so (|33|) appears to oscillate faster and faster as e — > 0, without a well-defined limit. 
Nonexistence of a pseudoclassical limit for the quantum dynamics was established in 36], by a stationary phase 
approach, with no recourse to the band formalism. It was nonetheless pointed out that, despite absence of such a 
limit, QAMs may be associated to certain rays, which do correspond to trajectories of some formally classical maps. 
The meaning of the latter maps is, at most, that of providing local phase space descriptions, near QAMs. Similar 
remarks apply in the case of the pseudoclassical maps (|27p . 

It is worth recalling that maps (|27p were derived from an ansatz, which would be optimally justified if the effective 
Hamiltonian of kicked dynamics could be replaced by the sum of the free and of the kicking Hamiltonians (sect. Illljl . 
This approximation is obviously invalid in a global sense, yet, in "spinless" cases, it is known to work remarkably well 
near stable fixed points 0] ; indeed, in the KR case it yields a pendulum Hamiltonian, which provides a good description 
of the motion near the stable fixed point of the Standard Map. This may be seen as a qualitative justification for the 
use of maps (|27[) . if restricted to the search of QAMs. 

D. Case q = 2. 

While the expressions in subsect. IIV Al are quite general, we may accomplish a detailed analysis when q — 2 and 
V(9) = kcos(8). In such a case the eigenvalues cOji'd, k) (j = 0, 1) of the resonant Hamiltonian can be written down 
explicitly (see appendix iBj) : 



— h (— 1) J arccos 



cos(/ccos(^/2) 



.4 v ' \ V2 

with m p — (—l) 2 ^ - . Therefore, our theory produces two maps (|28|) . which take the form: 



(34) 



J?i+1 =$t + Jt mod (2tt), 

' ^)) 

/l+sin2(/ccos(^±i)) 



/ „<7 \ sin ( k cos( , , 

J t+1 = J t +47rO-2(-l) J TO p fcsin(%M v L=LJd mod (2tt). (35) 



Going back to the time-dependent form the maps are written as 



d t+1 = -d t + AI t + 47r£li + g mod (27r), 

I t+1 =I t - {-iym p \ sin ^nZpL mod (2tt), 



(36) 

with q = 2 (eSj^i + nfl + t{3 — Trpfio) and where we have used ctj — —hfij,! (see appendix IB ip . 




FIG. 5: Probability inside a box of extension equal to L — 6 ~ AJq/\e\ in momentum, moving according to (I30|l for 
p — 3, q = 2, (r, s) = (1, 1) and e = —0.2828, as a function of quasimomentum /3 of the /?-rotor. A J is the size of the island in J. 
Dashed vertical lines refer to special values of quasimomenta, given by formula (|37p with f3o = v/3, v — 0, 1, 2, iVo = 0, m = 0. 
The probability is shown at time f = 100 (red) and t = 200 (blue). The parameter values are the same as in figlU 



We remark that in the case q = 2, avoided crossing between the eigenvalues (|34[) are absent for an arbitrary value 

Of*!. 

We may also check the selection criterion for quasimomenta, that in the present case assumes the form: 

Pi, v = --(N + fa) + J ° + - g + ft mod(l), (37) 
t qr A 

with /3q given by (Hi) in sect. [TTI and v = 0, ...,p — 1. A scan over possible /3 values reveals that indeed QAMs are 
greatly enhanced around the values predicted by (|32p : this is confirmed by fig0 in which the momentum probability 
transferred to the mode is shown as a function of 0. 



V. MODE SPECTROSCOPY AND CONNECTIONS WITH COLD ATOM EXPERIMENTS. 



A. Farey ordering of QAMs near a fixed resonance. 

We now elucidate how our findings apply to inspection of density plots like the one illustrated in figdl We point out 
that such a picture is of direct physical significance, since typical experimental protocols maintain k and g fixed, while 
performing a scan on the pulse period r. Such a scan, in the present context, has to be carried out around a resonant 
value, namely r(e) = 2-Kp/q + e. Density plots of momentum distribution disclose the presence of QAMs, since after a 
fixed number of kicks their momentum is linearly related to the acceleration (|30p : a depends on e, through the "bare" 
winding number qfl 

q n(e) = ^-g (2n P - + e) (38) 

and on the "dressed" winding number of the pseudo classical map r/s, which individuates the mode. We denote by 
n* the resonant (e = 0) value of the "bare" winding number (notice that for e = the maps correspond to pure 
rotation in J): 

il* = gfi(0) = 2tt^- 5 , (39) 
which is independent of the mapping index j. Formula (|39[) is a generalization of the analogous result found for q = 1 

USUI. 

As analyzed in [29T [52j for principal resonances, the parameter space of map (|28|) is characterized by the presence 
of regions (Arnol'd tongues), in which stable periodic orbits exist. Close to resonances we expect that mode- locking 
structure of the pseudo classical maps singles out modes whose winding number provide rational approximants to Q* : 
at the same time fat tongues are associated to small s values, so the corresponding modes should be more clearly 
detectable. This is the physical motivation underlying Farey organization of observed modes: whenever we observe 
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FIG. 6: Enlargement of fig[T]in the region 1.49 < t/2tt < 1.50625, around the resonance r res = 3n (p/q = 3/2). Momentum 
distributions is calculated after t — 200 kicks. Full lines show the theoretical curves (|30|) : the yellow ones refer to principal 
convergents of fl* . Starting from the left, the modes correspond to the stable periodic orbits of maps (|35p with: (r,s) = 
(14, 13), (25, 23), (12, 11), (23, 21) and (11, 10). 

two modes labelled by winding numbers ri/s% and r%/ ' S2 (r%/si < Q* < r-ij '32), the fraction with smallest denominator 
bracketed by the winding pair is the Farey mediant (r\ + r2)/(s\ + S2). 

We can now analyze in more detail figd which represents a numerical simulation of experimental momentum 
distribution after t = 100 vs r, which assumes values around a second order resonance, namely T rcs /2n = p/q = 3/2. 
All parameters are chosen to be accessible to experiments and the initial atomic distribution reproduces that employed 
in [22L [231. Izil. I25I] : g = 0.0386, k — 1 and the initial state is a mixture of 100 plane waves sampled from a gaussian 
distribution of momenta with FWHM ~ 9. 

Full lines in the figure delineate momentum profiles consistent with the acceleration (|30[) . with (r, s) given by 
winding number r/s of corresponding stable periodic orbits of maps (|28p . 

The value of Q* and the first few rational approximants (obtained upon successive truncation of the continued 
fraction expansion), corresponding to detectable modes, are: 



The first one (r, s) = (1, 1) is shown with yellow full line on the left of figffland the stability island of the corresponding 
fixed point is shown in fig[3] (a). The second and third are marked by full yellow lines in figEl which is an enlargement 
of fig[T]in the region 1.49 < r/27r < 1.50625, calculated for time t — 200. Farey organization is exemplified by the 
appearance of the (23,21) QAM, whose winding number is the Farey composition of the (11,10) and the (12,11) 
modes; the correspondent stable periodic orbit is plotted in figEtb). Through Farey composition law we may also 
identify observed modes to the right of r res , as shown in figEl 



The complexity of mode spectroscopy is further enhanced by the fact that, within some interval in r, arbitrarily 
many different resonant values occur. As a matter of fact it is possible to recognize in figffl modes coming from a wide 
set of resonances: besides q = 2 also q — 7, 15, 17, 21, 36, 40 contribute QAMs in the selected range; this is shown in 
fig[T]for q = 7 and in figJTJfor the other resonances. No QAM with q = 13 could be resolved in the range of figd] 

Farey composition is still of some use in the identification of the resonances to which modes belong: for instance, 
the very large mode on the right of the figure belongs to a QR between p/q = 3/2 and p/q = 2/1; applying Farey 
composition successively, we get the sequence p/q = 5/3, 8/5 (outside the plotted range in r) and then 11/7, to which 
the mode belongs. The accumulation point of the resonance p/q = 11/7 is ft* ~ 4.1923207 = 4 + [5, ...]. The mode 
shown in fig [T] corresponds to the first principal convergent of i.e. to the fixed point (r, s) — (4,1), shown in 
figEtc). The same occurs for the modes near resonances of higher q, shown in figEl 

We remark that a hierarchy in resonant fractions looks more cumbersome than the one considered for winding 
numbers, as for instance there does not seem to be any straightforward dependence on the size of q. Numerical 
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B. Visibility of resonances of different order. 
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FIG. 7: Enlargement of fig[T]in the region 1.45 < t/2tt < 1.49. The momentum distributions is calculated after t — 200 kicks. 
Full lines show the theoretical curves (|30|) : each color refer to a different quantum resonance, namely different values of p/q. 
Starting from the left the modes correspond to: p/q — 3/2 (theoretical curve not shown), 31/21,59/40,28/19,53/36,25/17 
and 22/15. These resonances p„/q„ lead to almost integer Q* and they are exacted from the sequence of Farey fractions 
obtained starting from po/qo = 1/1 and pi/qi = 3/2: ps/qs = 22/15, fig ~ 7.82566 = 8 — [5, 1, 2, 1, ...], (r, s) = (8, 1) (shown 
in black); p 9 /q g = 25/17, fig ~ 8.9165791 = 9 - [11,1,78,...], (r,s) = (9,1) (in purple); p 10 /q 10 = 28/19; fit ~ 10.0075 = 
10 + [131, ...], (r, s) = (10, 1) (in red); pu/qu = 31/21, 0*1 ~ 11.098678 = 11 + [10, 7, 2, 6, 1, ...], (r, s) = (11, 1) (in yellow). 
Further modes are shown in between the mentioned ones: p/q = 53/36 = 25/17©28/19 with Q* ~ 18.924152 = 19- [13, 5, 2, ...] 
(shown in pink) and p/q = 59/40 = 28/19 © 31/21 with tt* ~ 21.106256 = 21 + [9, 2, 2, 3, ...] (shown in orange). The (8, 1)- 
periodic orbit of the resonance ps/qs = 22/15 is plotted in figl^d). 



data however suggest that detectable modes appears in the vicinity of resonances leading to almost integer f2*, i.e. 
when the fractional part of O* is closer to the integers or 1 than to their Farey mediant 1/2. In these cases, the 
resonance may display a mode corresponding to a periodic orbit of period 1. As shown in fig[3 this condition may 
be fulfilled for different p/q values. Moreover, the absence of observable QAMs with q = 13 in the range of figlU 
even if p/q = 20/13 resonance belongs to the plotted t range, is consistent with this rough "thumb" rule. Indeed 
ft* ~ 7.4624908 = 7+ [2, 6, 6, ..], thus its fractional part is closer to 1/2 than to 0. 



The quantum dynamics of quantum accelerator modes, experimentally observed by exposing cold atoms to periodic 
kicks in the direction of the gravitational field, is theoretically described in terms of spinors, when the pulse period is 
close to a rational multiple of a characteristic time of the atoms (Talbot time). The reference model is a non-trivial 
variant of the well-known Kicked Rotator in an almost-resonant regime. If the detuning of the kicking period to 
the resonant value is assigned the role of the Planck constant, the problem is shown to share similarities with the 
semiclassical limit of the particle dynamics in presence of spin-orbital coupling. The separation of the spinor and 
orbital degrees of freedom, is based on an "adiabatic" assumption of Born-Oppenheimer type, valid for small detunings 
and for values of the parameters in which the QAMs manifest. In these parameter regimes, a description of some 
properties of the "slow" orbital motion, by means of formally classical equations, is finally achieved. Some results of 
a previously formulated "pseudo-classical" theory [HI, restricted to QAMs near principal resonances, are extended to 
arbitrary higher order resonances. Potential applications to current experiments on cold atomic gases are proposed. 

L.R. acknowledges useful discussions with Shmuel Fishman. 

APPENDIX A: SPINOR DYNAMICS AT EXACT RESONANCE AND WITHOUT GRAVITY. 

The Floquet operator at exact resonance in absence of gravity 77 = is given by (cp. Q) : 



VI. SUMMARY. 




- i7T E(N+/3 f 



(Al) 



where /3q = v/rp + rq/2 with c£Z. The evolution of spinor components under U Ies is given by |47j |: 

9-1 

4>M = 5> _ttt Vl e-^tC+ft)"^,^) (A2) 
(e-^) j _ i = (j| e -^«F e -^WFt e ^ § |/> = e "' (J "' )7 £ e -ia fm y_o e -ifc^(f+^ TO ) > (A3) 

" m=0 

F is the Fourier transform in C 9 : (i|F|Z) = e^ l ~^ 1 / \fq and S is the spin operator Hence from (|A2|) . the resonant 
Floquet operator is decomposable in spin propagators A(i?), given by unitary matrices of rank g; namely it acts as 

(£W)(0) = A(0#(i?). 

APPENDIX B: CASE g = 2: RESONANT EIGENVALUES AND EIGENVECTORS. 

At primary 2nd order (q = 2) resonances, the resonant values of quasimomentum are /3q = v/p with v — 0, 1, ..,p—l. 
We choose /?o = 0. From ((HJ) and denoting w(t9; fc) = fccos($/2) and m p = (— l)^ - , we find 

A k 20= f cos(v($;k)) m p e 4 f sm(v(d;k)) \ ^ 
' p ' ' \ — je~ 4 2 sin(u('i3; k)) im p cos(v (i?; fc)) / 

Matrix (jBip can be written in terms of Pauli matrices as follows: 

A fe , p , 2:0 = e im ^ (cos(w(0;fc))+isc(0;A;)-*), (B2) 
= e lm ^ (cos(tf(i?;fr))+is^w(i?;fc)) (n(?9; fc) • <x)) , (B3) 

(B4) 



£j(tf; fej = arccos ^ -j= I 



where the vector cc(i?; fc) E M 3 has components: 



a;i(i?; fc) = m p sin I — — Wp-j J sin(t>($; fc)), 
x 2 (i}; k) = sin I - + ro p - J sin(u(tf; fc)), 

77? 

x 3 (tf;/c) = ^cos(v(??;fc)); (B5) 

V 2 



n = i/a; is a unit vector in K 3 , and cc(i?; fc) = | |jc(i?; 

Using a well-known formula, equation (|B2[) may be written: 



A fc ,p, 2 ,o = e im "^e iQ ^ (B6) 
which directly yields to the Resonant Hamitonian for this case: 

HS, 2 ,o = - [m p J + fe) fc) • <r] . (B7) 
As matrix n ■ <x has eigenvalues ±1, the eigenvalues of the resonant fiber are: 

A w (tf,fc) = e ™ P [f+(-i)M<?;fc)] = e -«^(M) ( B8 ) 

cjjft k) = -mp [J + (-l)^(^; k)] j = 0, 1. (B9) 



Normalized eigenvectors are: 



^^e^O^f; (BIO) 

yi M)^(Vf7 (Bll) 

aj(tf;k) = 1 + sin 2 (v(i?;fc)) + (-1) J ' cos(v(tf; fc))^/l + sin 2 (z;(i9; fc)), (B12) 

fc) = cos(w(i9; fc)) + (-1)^1 + sin 2 (u(z9; fc)). (B13) 

The 7j(i?) are arbitrary phases. For uqtt < k < (no + 1)tt (no 6 N), ai($;fc) has a finite set of zeros in {$„ = 
±2 arccos(-7m/fc), < n < no}. Discontinuities in the zeros of ai(#; k) may be removed by appropriate choices of 71. 
For k < 7r and — 7r < ?9 < 7r, one may choose, e.g. 7j(#) = f <5j,i- 

For increasing values of the kicking strength fc, eigenvalues (jB9|) display thicker and thicker oscillations. An example 
is shown in figJ2](a), for k = 1,3,5. Nevertheless, contrary to higher values of q, in case q = 2, the two eigenvalues 
(|B9j) neither cross nor become closer that a minimal gap, equal to 7r/2, for arbitrary high values of fc. As a matter of 
fact, the band width of each eigenvalue, defined as Bj(k) = \u>(-d m ax', k) — £>(i? m j n ; k)\, with z9 ma2; and $ m in absolute 
maximum and minimum points in [0, 2w[, does not exceed 7r/2. For k < ir, Bj(k) is an increasing function of equal 
to Bj(k) = \u(0]k)-u>(ir,k)\ = |arccos(cos(fccos(?9/2)/v / 2)) -tt/4|; for k > it, Bj(k) = |<D(i?i; fc) - u>(tt; fc)| = vr/2. 

For /3q = ^/p, the same results hold, apart for a constant phase factor in (|Blj) and m p replaced by — m p for i/ odd. 



1. Vector and Scalar potential. 



The vector and scalar potentials can be explicitly computed, using analytical expressions of the eigenvectors of 
resonant fiber (and choosing phases 7j(i?) = § £7,1 )■ First one computes: 



1 



2 a 



4w = E z2 i(^iW)i 2 = ^'i' 



1=0 
1 



;=n 3 \ j / 



»<^i(i?)|^W> = "7iW 



1%! 
4 oj 



(^•wi^W) = - (aYcoJv + hbtf + l( b rf - iM) + \ Sjil (1- M 



Oj V 2 8 4 «J / 4 



and then the vector and scalar potentials are given by: 



^iW = ai=-^i,i, (B14) 



B M = _L ( {i f cos * v + i(6,) 2 - 7%^ = „ (B15) 



4 Oj / (l + sin 2 i;) 2 ' 

with v ee fc) = k cos(i?/2). 

For g = 1 the resonant fiber is a matrix with a single element exp(— iv($, k)), with eigenfunction </?o($) 
exp(ini?)/v27T (n € Z); 6*00 = and therefore (f2"Tj) it yields -4q(i9) = ao = — ti with n € Z and Bo(^) = 0. 



[1] G. Casati, B. V. Chirikov, F. M. Izrailev, and J. Ford, in Lectures Notes in Physics , edited by G. Casati and J. Ford 
(Springer, Berlin, 1979), Vol. 93. 



[2] F.M.Izrailev, Phys. Rep. 196, 299, (1990); S. Fishman, in Proceedings of the Inter-national School of Physics Enrico Fermi; 

Varenna Course CXIX, G.Casati, I.Guarneri and U. Smilansky eds. (North- Holland 1993). 
[3] R. Lima and D. Shepelyansky, Phys. Rev. Lett. 67, 1377, (1991). 
[4] T. Geisel, R. Ketzmerick, and G. Petschel, Phys. Rev. Lett. 67, 3635, (1991). 

[5] R. Artuso, F. Borgonovi, I. Guarneri, L. Rebuzzini and G. Casati, Phys. Rev. Lett. 69, 3302, (1992). 

[6] I. Guarneri and F. Borgonovi, J.Phys. A Math. Gen. 26, 119, (1993). 

[7] A. Backer, R.Ketzmerick, and A. Monastra, Phys. Rev. Lett. 94, 054102, (2005). 

[8] A. Backer, R.Ketzmerick, S. Lock and L. Schilling, Phys. Rev. Lett. 100, 104101, (2008). 

[9] M. Sheinman, S. Fishman, I. Guarneri, L. Rebuzzini, Phys. Rev. A 73, 052110, (2006). 
[10] J. Feist, A. Backer, R. Ketzmerick, S. Rotter, B. Huckestein, and J. Burgdo, Phys. Rev. Lett. 97, 116804, (2006). 
[11] A. Ossipov, T. Kottos, and T. Geisel, Europhys. Lett. 62 (5), 719, (2003). 
[12] F. Borgonovi, I. Guarneri, and L. Rebuzzini, Phys. Rev. Lett. 72, 1463, (1994). 

[13] F.B. Dunning, J. C. Lancaster, C. O. Reinhold, S. Yoshida and J. Burgdorfer, Adv. At. Mol. Opt. Phys 52, 49, (2005). 

[14] G.G. Carlo, G. Benenti and D.L. Shepelyansky, Phys. Rev. Lett. 95, 164101, (2005). 

[15] H. Schanz, M.F. Otto, R. Ketzmerick, and T. Dittrich, Phys. Rev. Lett. 87, 070601, (2001). 

[16] C. E. Creffield, G. Hur, and T. S. Monteiro, Phys. Rev. Lett. 96, 024103, (2006). 

[17] I. Dana, V. Ramareddy, I. Talukdar,and G. S. Summy, Phys. Rev. Lett. 100, 024103, (2008). 

[18] F.L. Moore, J.C. Robinson, C.F. Bharucha, B. Sundaram, and M.G. Raizen, Phys.Rev.Lett. 75, 4598, (1995). 

[19] H. Amman, R. Gray, I. Shvarchuck, N. Christensen, Phys. Rev. Lett. 80, 4111, (1998), 

[20] P. Szriftgiser, J. Ringot, D. Delande, J.C. Garreau, Phys. Rev. Lett. 89, 224101, (2002). 

[21] C. Ryu, M.F. Andersen, A. Vaziri, M.B. d'Arcy, J.M. Grossmann, K. Helmerson, W.D. Phillips, Phys. Rev. Lett. 96, 
1604031, (2006). 

[22] M.K. Oberthaler, R.M. Godun, M.B. d'Arcy, G.S. Summy and K. Burnett, Phys.Rev.Lett. 83, 4447, (1999). 

[23] R.M. Godun, M.B. d'Arcy, M.K. Oberthaler, G.S. Summy, and K. Burnett, Phys. Rev. A 62, 013411, (2000). 

[24] S. Schlunk, M.B. d'Arcy, S.A. Gardiner, G.S. Summy, Phys. Rev. Lett. 90, 124102, (2003). 

[25] M. B. d'Arcy, R.M. Godun, D. Cassettari, and G. S. Summy, Phys. Rev. A 67, 023605, (2003). 

[26] M.V.Berry and E.Bodenschatz, J. Mod. Opt. 46, 349, (1999). 

[27] S. Fishman, I. Guarneri, and L.Rebuzzini, Phys. Rev. Lett. 89, 084101, (2002). 

[28] S. Fishman, I. Guarneri, and L.Rebuzzini, J. Stat. Phys. 110, 911, (2003). 

[29] I. Guarneri, S. Fishman, and L. Rebuzzini, Nonlinearity 19, 1141, (2006). 

[30] R. Hihinashvili, T. Oliker, Y.S. Avizrats, A. Jomin, S. Fishman, and I. Guarneri , Physica D 226, 1, (2007). 
[31] G. Behinaein, V. Ramareddy, P. Ahmadi, and G. S. Summy, Phys. Rev. Lett 97, 244101, (2006). 
[32] L. Rebuzzini, R. Artuso, S. Fishman, and Italo Guarneri, Phys. Rev. A 7 6, 031603(R), (2007). 
[33] P. L. Halkyard, M. Saunder s S. A. Gardiner, and K. J. Challis, \arXw:0807.2587\ (2008). 

[34] G. Lemarie and K. Burnett, {arXiv:quant-ph/0602204\ (2006). 

[35] M. Saunders, P. L. Halkyard, S. A. Gardiner and K. J. Challis, \arXiv:0806.389j\ (2008). 
[36] I.Guarneri and L. Rebuzzini, Phys. Rev. Lett. 100, 234103, (2008). 
[37] G. Summy, private communication. 

[38] F.M. Izrailev and D.L. Shepelyansky, Theor. Mat. Phys. 43, 353, (1980). 
[39] R.G. Littlejohn and W.G. Flynn, Phys. Rev. A 44, 5239, (1991). 
[40] R.G. Littlejohn and W.G. Flynn, Phys. Rev. A 45, 7697, (1992). 
[41] G. Casati and I. Guarneri, Comm. Math. Phys. 95 , 121, (1984). 

[42] V.V. Sokolov, O.V. Zhirov, D. Alonso, and G. Casati, Phys. Rev. Lett. 84, 3566, (2000). 
[43] V.V. Sokolov, O.V. Zhirov, D. Alonso, and G. Casati, Phys. Rev. E 61, 5057, (2000). 
[44] J. Farey, Philos. Mag. 47, 385, (1816). 

[45] G.H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, (Clarendon, Oxford, 1979). 
[46] I. Dana and D.L. Dorofeev, Phys. Rev. E, 73, 026206, (2006). 
[47] S.J. Chang and K.J Shi, Phys. Rev. A 34, 7, (1986). 

[48] S. Teufel, Adiabatic Perturbation Theory in Quantum Dynamics, Lecture Notes in Mathematics 1821, Springer 

Berlin/Heidelberg, (2003). 
[49] I. Dana, E. Eisenberg, and H. Shnerb, Phys. Rev. Lett. 74, 686, (1995). 
[50] M.V. Berry, Proc. R. Soc. A 392, 45, (1984). 
[51] B. Simon, Phys. Rev. Lett. 51, 2167, (1983). 

[52] A. Buchleitner, M.B. d'Arcy, S. Fishman, S.A. Gardiner, I. Guarneri, Z.Y. Ma, L. Rebuzzini, G.S. Summy, Phys. Rev. 
Lett. 96, 164101, (2006). 



